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Abstract 

A preliminary measurement of \V c b\ and the branching fraction B(B° — * D* + l~V() has been per- 
formed based on a sample of about 55,700 B° — ► D* + i?~z^ decays recorded with the BABAR detector. 
The decays are identified in the D* + — > D°tt + final state, with the D° reconstructed in three dif- 
ferent decay modes. The differential decay rate is measured as a function of the relativistic boost 
of the D* + in the B° rest frame. The value of the differential decay rate at 'zero recoil', namely 
the point at which the D* + is at rest in the B° frame, is predicted in Heavy Quark Effective 
Theory as a kinematic factor times .F(l)|Vd>|, where T is the unique form factor governing the 
decay. We extrapolate the measured differential decay rate to the zero recoil point and obtain 
T(\)\V ch \ = (34.03 ± 0.24 ± 1.31) x 10~ 3 . Using a theoretical calculation for .F(l) we extract 

|V c6 | = (37.27 ± 0.26(stat.) ± 1.43(syst.)li j(theo.)) x 10~ 3 . 

From the integrated decay rate we obtain 

B(B° -» D* + rv£) = (4.68 ± 0.03 ± 0.29)%. 
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1 Introduction 



In the framework of the Standard Model, the weak coupling between quarks of different flavors is 
described by the Cabibbo-Kobayashi-Maskawa (CKM) matrix. Its elements are not predicted by 
theory, but are constrained only to give a unitary matrix. 

A precise measurement of \V c b\, the element corresponding to the b —> c transitions, is needed 
to determine whether the CKM matrix provides a quantitatively accurate description of the CP 
violation observed in B° mesons. Progress in the phenomenological description of heavy flavor 
semileptonic decay allows the determination of \V c b\ with small theoretical uncertainty, either from 
the inclusive process b — > civ, or from an analysis of the form factors in the decay B° — > D* + £~i7£. 
The present measurement is based on the second approach [1]. 

The decay rate for the process is proportional to | V& | 2 and to the square of the hadronic matrix 
elements describing the transition from a B° to a D* + meson. In the context of the Heavy Quark 
Effective Theory (HQET) [2], the matrix element is proportional to a single form factor T(w), 
where w is the product of the B° and D* + four- vector velocities. For w = 1, the D* + is produced 
at rest in the B° rest frame. Heavy-flavor symmetry predicts the normalization .F(l) = 1 in the 
limit of an infinitely massive 6-quark. Corrections to this prediction due to perturbative QCD have 
been computed up to second order in a s [3] and the effect of finite b and c quark masses has been 
calculated in the framework of HQET, yielding = 0.913^;^ [4]. 

A measurement of the differential decay rate near w = 1 determines \ V c b\ with small theoretical 
uncertainty. The rate at w = 1 is suppressed by phase space. Consequently, the differential rate 
dT/dw is measured, where F(w) is parametrized using several different functional forms [5, 6] (see 
also discussion below), and is extrapolated to w = 1. Results based on this approach have been 
reported by the ARGUS [7], Belle [8], and CLEO [9] collaborations operating at the T(4S) and by 
ALEPH [10, 11], DELPHI [12, 13] and OPAL [14] at LEP. 

This paper presents preliminary results for the measurements of \V c b\ and of the branching 
fraction B(B° — > D* + l~V() performed on a sample of B° — > D* + £~t7£ decays selected from events 
with a high momentum charged lepton i~ and a D* + 1 . The D* + is reconstructed from its decay to 
a charged pion and a D°. This pion is produced with small momentum, and is commonly referred 
to as the slow pion (vr s i ow ). 

This paper is organized as follows. The next section is dedicated to a brief description of the 
BABAR detector and of the data sets employed. Section 3 describes the event selection and the 
composition of the sample used for the measurement. The fit method and its results are described 
in Section 4; the evaluation of the systematic error is discussed in Section 4.4. Finally, Section 5 
presents the conclusions and the comparison with other measurements of |V^,|. 

2 The BABAR Detector 

The data were collected with the BABAR detector at the PEP-II asymmetric energy e + e~ storage 
ring [15]. The sample consists of 79.1 fb™ 1 collected at a center-of-mass energy corresponding to 
the T(4S) resonance ( "on-resonance" ), and 9.6 fb" 1 collected 40 MeV below the T(4S) resonance 
("off-resonance") for continuum background studies. The on-resonance sample corresponds to 
N T = (85.9 ± 0.9) x 10 6 T(4S) decays to BB mesons. 

The BABAR detector is described in detail elsewhere [16]. The momenta of charged particles 
are measured using a combination of a five-layer silicon vertex tracker (SVT) and a 40-layer drift 

1 Charge conjugate states are always implicitly considered. Lepton as used here means either electron or muon. 
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chamber (DCH) in a 1.5 T solenoidal magnetic field. A detector of internally-reflected Cherenkov 
radiation (DIRC) is used for charged-particle identification. Kaons are identified with a neural 
network based on the likelihood ratios calculated from dE/dx measurements in the SVT and DCH, 
and from the information from the DIRC. A finely segmented CsI(Tl) electromagnetic calorimeter 
(EMC) is used to detect photons and neutral hadrons, and to identify electrons. Electron candi- 
dates are required to have a ratio of EMC energy to track momentum, EMC cluster shape, DCH 
dE/dx, and DIRC Cherenkov angle all consistent with the electron hypothesis. The instrumented 
flux return (IFR) contains resistive plate chambers for muon and long-lived neutral hadron identifi- 
cation. Muon candidates are required to have IFR hits located along the extrapolated DCH track, 
and energy deposition in the EMC and penetration length in the IFR consistent with a minimum 
ionizing particle. 

3 Event Selection and Sample Composition 

3.1 Event Selection 

Events are selected that contain a fully-reconstructed D* + and an oppositely charged lepton. Those 
D* + (r combinations not coming from the signal decay are suppressed by means of kinematic cuts. 
A candidate lepton is searched for among all charged tracks with momentum greater than 

1.2 GeV/c in the T(4S) rest frame. Electron candidates are selected with an efficiency of about 
90% and a hadron misidentification probability of less than 0.2%. Muon candidates are selected 
with an efficiency of ~ 60% and hadron misidentification probability is about 2%. 

The D* + candidate is selected in the decay mode D* + — > D°ir + and the D° meson is recon- 
structed in the three modes K~Tr + , K~ir + ir + ir~ , and K~ir + 7r°. The tt° is reconstructed from two 
photons, each with energy larger than 30 MeV, and must have a total energy larger than 200 MeV 
and an invariant mass between 119.2 and 150.0 MeV/c 2 . The invariant mass of the two photons is 
constrained to the 7r° mass and the pair is kept as a 7r° candidate if the x 2 probability of the fit 
exceeds 1%. Charged kaon candidates satisfy loose identification criteria for the K~ir + mode and 
tighter criteria for the K~7r + 7r + ir~ and K^tt^tt modes. For the K~it + it mode, the D° candidate 
is retained if the square of the decay amplitude in the Dalitz plot for the three-body candidate, 
based on measured amplitudes and phases [17], is larger than 10% of its maximum value. D° 
candidates are accepted if they have an invariant mass within 17 MeV/c 2 of the D° mass for the 
K~it + and K~it + it + it~ modes, and within 34 MeV/c 2 for the K~it + it° mode. The invariant mass 
of the decay products is then constrained to the D° mass [18] and the tracks are constrained 
to a common vertex by a simultaneous fit. The D° candidate is retained if the x 2 probability of 
the fit is more than 0.1%. The low-momentum pion candidates for the D* + — > D vr^ decay are 
selected from among tracks with total momentum in the laboratory frame less than 450 MeV/c, 
and momentum transverse to the beam line in the laboratory frame greater than 50 MeV/c. The 
momentum of the D* + candidate in the T(4S) rest frame is required to be between 0.5 and 2.5 
GeV/c. These requirements retain essentially all signal events and reject higher momentum D* + 
from continuum events. Candidates from D* + are preselected with the cut on the mass difference 
AM = M D o,r slow - M D o < 165 MeV/c 2 . The AM distribution has a kinematic threshold at the 
mass of the 7r + , and a peak at 145.5 MeV/c 2 with a resolution of about 1 MeV/c 2 . The combina- 
toric background below the signal is evaluated using events in the AM sideband. A harder AM 
cut is applied later. 

The x 2 probability of the fit to a common vertex of the lepton, the 7r s i ow , and the D° candidate, 
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constrained to the beam spot, must be greater than 1%. 

Candidates for D* + l~ are selected if | cos A#* hms J < 0.85, where A#* hrust is the angle between 
the thrust axis of the D* + l~ candidate and the thrust axis of the remaining charged and neutral 
particles in the event. The distribution of | cos A#^ irust | is peaked at 1 for jet-like continuum events, 
and is flat for BB events, where the two B mesons decay independently and practically at rest with 
respect to each other. This cut retains about 85% of signal and rejects about 47% of continuum 
events. 

Finally, the angle between the lepton and the D* + in the T(4S) rest frame must satisfy 
cos0£,*£ < 0. This cut reduces significantly the events due to random combinations of a lepton and 
a D* + produced by two different B mesons, while keeping most of the signal events. 

3.2 Sample Composition 

Several processes contribute to the selected sample: 

1. Signal decays. 

2. Combinatoric background from any random combination of tracks mimicking a true D° or 
from a correctly identified D°, combined with a low- momentum charged track not originating 
from D* + decay. This category includes events from BB and from continuum production. 
Because of its combinatorial origin, the AM distribution does not exhibit a peak. All the 
other backgrounds listed below peak in AM . 

3. Fake-lepton background, in which a charged hadron is misidentified as a lepton and is com- 
bined with a correctly reconstructed D* + . Due to the very good performance of the electron 
identification, the rate of fake electrons is negligible. 

4. Continuum background, where a true D* + and a lepton (true or fake) are produced from 
e + e~ — ► cc events. 

5. "Uncorrelated" background, in which the lepton and the D* + originate from two different 
B mesons. Most of these combinations are removed by the requirement that the lepton and 
D* + have opposite charge. Opposite charge combinations arise mainly from events with B°B° 
oscillations where both B mesons decay into the same flavor eigenstate, or events in which 
the D* + comes from the hadronization of the virtual W from B decay, and thus has a charge 
opposite to the one expected from its parent B. 

6. B — ► D*^ 't~ V(_X background, where one (or more) hadrons are produced in addition to the 
D* + , for either neutral or charged B. A considerable fraction of this background is due to 
the intermediate production of orbitally excited D states, which then decay to a D* + ir final 
state; this is referred to as D** background in the following. 

7. "Correlated" background, in which the B decays to a D* + and a heavy particle (either a charm 
meson or a r), which then decays to a lepton. These events have the same charge correlation 
as signal events, but are suppressed by their low branching fraction and by the kinematic 
cuts. Such events correspond to a few percent of the peaking sample. Their contribution is 
fixed to the value computed in the simulation. 



10 




140 145 150 155 160 165 140 145 150 155 160 165 

D*-D° Mass Difference (MeV) D*-D° Mass Difference (MeV) 



Figure 1: AM distribution for events passing all selection criteria: D* + e~ (left) and D* + fi~ 
(right) candidates in on-resonance data. Points correspond to data, and the curve is the result of 
unbinned maximum likelihood fits to this signal sample, the fake-lepton control sample and the 
rescaled off-resonance events. The shaded distributions correspond to the continuum, fake-lepton, 
and combinatorial backgrounds described in the text. The correlated-lepton, the B — > D* + d.~i7^X 
and the uncorrelated backgrounds are included in the signal distribution. 



The fraction of combinatoric events below the D* + peak is determined from a fit to the AM 
distribution in the region 0.138 < AM < 0.165 GeV/c 2 (see Figure 1). In the fit, D* + events are 
described by the sum of two Gaussian functions while the shape of combinatoric events is described 



by the function 1 — exp ( — Am m7r ) (£2) , where the overall normalization, c\, and C2 are 



Cl 



Am 



(-'_> 



determined by the fit. The width of the AM distribution for signal events is dominated by the 
experimental resolution. In one third of events the vr s i ow track is reconstructed in the SVT and in 
the DCH; in the remainder, the track is reconstructed only in the SVT. As the resolution in the 
first sample is better, the two samples are fitted separately. 

The AM distribution of a fake-lepton control sample is fitted simultaneously to determine the 
fake-lepton background. The control sample is obtained by selecting events with the kinematic cuts 
described above, except for the requirement that the candidate lepton fails a loose lepton selection 
criterion. Lepton identification efficiencies and hadron misidentification probabilities, estimated 
from data control samples of pure leptons and hadrons, are used to scale the number of peaking 
events in this hadron-substituted sample to the expected amount of fake-lepton background in the 
signal sample. 

The continuum contribution is fitted from the off-resonance event sample and scaled according 
to the ratio between the on-resonance and off-resonance luminosities. 

For each D° final state, the signal, fake-lepton and continuum samples are fitted simulta- 
neously. The mean values, widths and the relative normalization of the two Gaussian func- 
tions describing the signal are common to the three data sets, while the parameters describ- 
ing the shape of the combinatoric background are fitted independently to each data sample. 
Cuts on AM are then applied to reduce the amount of combinatorics in the subsequent stages 
of the analysis: 0.143 < AM < 0.148 GeV/c 2 if the 

"^siow is reconstructed by the SVT alone, 

0.144 < AM < 0.147 GeV/c 2 otherwise. 
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The fraction of uncorrelated and D** events is determined by exploiting the decay kinematics. 
The angle between the B° and the pseudo-particle obtained by adding the D* + and the t~ four- 
momenta is computed according to the following expression: 

a -(M| + - 2E B0 E DH ) + M^ iss 

cosO B o DH = — — , (1) 

*Pb°Pdh 

where M m i ss is the invariant mass of all the other particles produced in association with the lepton 
and the D* + in the decay of the B°. We compute cosd B o D *^ by fixing M miss to zero, under the 
hypothesis that this is a true B° — > D* + l~i/£ decay, for which only a neutrino is missing. Neglecting 
the broadening due to experimental resolution, the obvious constraint — 1 < cos#g<> < 1 applies 
to signal events. Due to the production of additional particles, events from D** background have 
positive M^ iss values and therefore produce large tails below -1 in cos 9 bo d*i distribution, while 
uncorrelated background events preferentially populate the region with cos 9 B o £,*£ > 1 . This 
allows for the determination of the fraction of D** , uncorrelated, and signal events in the sample, 
by means of a fit to the cos9 B o jB .*£ distribution in data (see for instance Figure 2). The shapes of all 
components in the fit are determined from Monte Carlo simulation, while the number of events for 
each of the other background categories is determined from the fit to the AM distributions. The 
same procedure is applied to 18 independent samples, identified according to lepton flavors (x2), 
period of data acquisition (x3), D decay mode (x3). The data sets are further split into bins 
of w (defined in Section 4.1) in order to reduce the systematic error corresponding to variations 
in background composition and normalization with w. The resulting large statistical uncertainties 
are reduced by requiring that the fractions of uncorrelated and D** background with respect to 
the signal be the same for electrons and muons, and for all D° decay modes. The validity of 
these hypotheses is verified by Monte Carlo simulation. The distributions in cos9 B o D *£ for the 
B° -► D* + e-i7 e ,D° -► K~tt + and B° -► D* + pT u^, D° -► K~ir + ir° samples (integrated over w ) 
are shown in Figure 2. 

Only events with | cos 9 B o D *£\ < 1.2 are further analyzed. This selects a final sample of 55700 
signal events for the measurement. 



4 Determination of \V c b\ 

4.1 Parametrization of the Decay Width 

Assuming the connection between the various form factors provided by HQET, the expected number 
of signal events can be expressed as a function of w by the relation 



cbV 
dw 

dB 

dw 



= 4 N T /oo B(D*+ - D\+ ow ) B(D° - K~mr) e(w) 



C 2 

' ! mL+(M b » - M D ,+ ) 2 Vw 2 - 1 (w + l) 2 



dB 

dw ' 



48vr 3 t B o " D 



Vd, | 2 F 2 (w) 



(2) 



1 + 



Aw 1 — 2wr + r 2 
TTw (l-r) 2 



The factor of 4 accounts for the fact that two B mesons are produced in each event, and that both 
electrons and muons are used; Nr is the number of T(4S) produced, /oo is the fraction T(T(AS) — > 
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cos ®BO,D*l 



Figure 2: Distributions in cosObo D*e f° r the B° —> D* + e v e ,D°^K ir + (left) and 
B° —*■ D* + pTu^ D° — ► K~ ~tt + it samples (right). Points with error bars show the data. Shaded 
areas represent the Monte Carlo distributions. The fit determines the normalization of the Monte 
Carlo distributions. The tail to the left in the signal electron sample is due to bremsstrahlung. 
Below the histograms are shown the fractional deviations of the data from the fit results. 

B°B°)/T(T(AS)) 2 , B(D* + D°ir+ l(m ) is the branching fraction for the decay D* + D°ir + , 
B(D° — > K~nir) is the branching fraction for the decay of the D° to the final state considered, t b o 
is the B° lifetime and r is the ratio of meson masses, r = M D *+/M B o. The values employed for 
these parameters, as determined by independent measurements, are reported in Table 1. 

Table 1: Values of inputs to the decay width calculation, as obtained from Ref.[18]. 



Parameter 


Value 


B{D* + -> D U 7T+) 
B(D° -> K-7T+) 
B(D° -> K-iT+iT+iT-) 
B(D° -> K-it + it ) 

f+-/f 00 


67.7 ± 0.5% 
3.80 ± 0.09% 
7.46 ± 0.31% 
13.1 ± 0.9% 
1.548 ± 0.018 ps 
1.055 ±0.055 



The iu-dependent reconstruction efficiency e(w) is determined by means of the detailed detector 
simulation. The w dependence of the form factor J~[w) is unknown. Since only a small range of 
w is allowed by phase space, a Taylor series expansion limited to second order has typically been 

2 The value /oo = 0.487 ± 0.013 is obtained from Table 1 assuming that the T(4S) decays to only to charged or 
neutral B mesons, i.e. imposing the constraint /oo + /h — = 1. 
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used: 



F{w) = (1 + P %{1 -w) + c(l - w) 2 + 0(1 - wf). (3) 

Apart from .F(l), the theory does not predict the values of the higher order coefficients, which must 
be determined experimentally. The first measurements of |V^,| were performed assuming a linear 
expansion, i.e., setting c = ([7, 9, 10, 13]). However, the requirement of analyticity and positivity 
of the QCD functions describing the local currents leads to the prediction that c should be positive 
and it should be related to the heavy meson radius Pjr{see [5]) by the relation: 



0.66/4- - °- n • 



(4) 



Results employing this analyticity bound were obtained by the ALEPH, DELPHI and OPAL 
collaborations (see [11, 12, 14]). An improved parametrization based on the above consideration is 
proposed in [6], accounting for higher order terms and so reducing to ±2% the relative uncertainty 
on 1 14ft | due to the form factor parametrization. A new function Ai(w) is introduced in that 
calculation, connected to F{w) by the following relation: 



F 2 (w) x 
A\{w) x 



1 + 



4u> 1 — 2wr + r 
1 + w (1 — r) 2 



(5) 



1 — 2wr + r 
(1 - r 



2 r 



1 + 



w 



1 



w + 1 



Ri{wf 



+ 



1 + 



w 



1 2 



R2H) 



where Ri(w) and R2(w) are ratios of axial and vector form factors also given in Ref. [6]. The 
following parametrization, depending only on a single unknown parameter p 2 ^, is proposed for 
Ai(w): 



Ai(w) = Ai(l) x 



8 / 9 2 4l z + (53p 2 4l 



I5)z 2 - (231/4 - 91)z 2 



with 



_ VwTT- V2 

\/w~+l + y/2' 

It should be noted that, in the limit w — > 1, Ai(w) — > J-(w), so that Ai(l) = ^"(1). 

Monte Carlo events employed for this analysis were produced with a linear parametrization for 
T{w), while experimental data are fitted using expressions of Ref. [6]. 

Experimentally the variable w can be expressed as: 



w 



(M| t ,+M 2 , + -g 2 ) 
(2M B oM D ,+ ) 



where q 2 = (p B o ~Pd*+) 2 . The magnitude of the B° momentum is known; its direction is obtained 
from Equation 1 with an azimuthal ambiguity about the direction of the D* + £~ pair. The two 
extreme solutions corresponding to the minimal and maximal angles between the B° and the D* + 
are used to define the quantity 



+ w n 



w 



The simulation shows that w is a good estimator of w, with a resolution a(w — w) 
sponding to about 4% of the full physical range. 



(6) 

0.02, corre- 
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4.2 Fit Method 



Data and Monte Carlo simulated events are collected in ten equal-size w bins. A few events, which, 
due to resolution, have w > w max = 1.503 are discarded. A least-squares fit is then performed 
comparing the number of events observed in each bin to the sum of signal and background events. 
The number of background events in each bin is determined as explained above and is fixed in the 
fit. The signal contribution is obtained at each step in the minimization procedure by properly 
weighting each generated event surviving the selection. The \ 2 i s: 

!\ri 

10 (Kfi _ i\ri _ v MC W\ 2 
*=! ^data + °bck + Ej=l C 

where the index i runs over the ten w bins; ^datafMC) is the number of data (signal Monte Carlo) 
events found in the i th bin; N^ ck , a^ ck are the numbers of background events and their errors. The 
weight for each Monte Carlo event is computed as the product of three terms 

Wj = W c x W] x W] heo , (8) 

where 



• W 1 -' is an overall fixed scale factor, that accounts for the relative luminosity of data and 
signal Monte Carlo events, for the difference in the actual branching fractions relevant to this 
analysis 3 , and any possible overall efficiency scale factor. 

• Wf accounts for the efficiency correction, computed event per event as the ratio of the ef- 
ficiency in the real data and in the Monte Carlo. It accounts for differences in tracking, 
particle identification, 7r° reconstruction, and depends on the three-momenta of the lepton, 
the soft pion and of all the particles from the D decay. It remains unchanged at all steps 
of the minimization. (If the Monte Carlo simulation were perfectly tuned, this factor would 
be identically one). It should be noted that only the factor due to the slow pion tracking 
efficiency introduces a net dependence of this term on w. 

• W\ heo = ftheo(w, p 2 , \ V c b\)/ /mc(w; py C , IV^Imc) is the term accounting for the theoretical 
function describing the decay. It depends on the parameters to be determined («4i(l)|Vcb| 
and p\^) and varies for each step of the minimization. 

Here the function /theo corresponds to the expressions in Equations 3 and 5, while the function 
used to generate Monte Carlo events, /mc, contains a simpler, linear parametrization of the form 
factor dependence on w. 

The value of the branching ratio is then computed by integrating the differential expression in 
Equation 3. 

The fit has been performed as a blind analysis, i.e., the values of „4i(l)|V c b| and were hidden 
until the study of the systematic errors was completed and all consistency checks were performed. 

3 obviously not including B(B° -> D* + £~i7 e ) 
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Table 2: Fit results and the branching ratio B(B° — > D* + l~V() , for sub-samples defined by D° 
decay mode, lepton species, and year of data taking. The third columns shows the statistical 
correlation between *4i(l)|Vd,| and p\ x . 





L> u /lepton 


Ai{l)\V ch \ x 10 3 


pX 


Govr stat 


B % 


nx z ) % 




Kir/e 


34.69±0.95 


1.284± 0.083 


91% 


4.72± 0.11 


3.1 




Kmnr/e 


35.68±1.42 


1.452± 0.106 


93% 


4.55± 0.16 


98.4 


o 
o 


Kmr°/e 


32.36±1.17 


1.241± 0.109 


93% 


4.21± 0.12 


20.7 


o 

CM 


Kir/fi 


33.00±0.99 


1.145± 0.097 


92% 


4.61± 0.11 


21.6 




Ktttttt/ fl 


33.51il.48 


1.410± 0.127 


93% 


4.11± 0.15 


26.0 




Kmr / fi 


32.75il.20 


1.235± 0.111 


93% 


4.32± 0.13 


83.1 




Kir/e 


34.92±0.73 


1.237± 0.063 


92% 


4.91± 0.08 


6.8 




Kmnr/e 


35.06±1.08 


1.357± 0.085 


93% 


4.63± 0.12 


45.6 


i—l 

o 


Kmr°/e 


33.68±0.89 


1.190± 0.081 


93% 


4.69± 0.09 


70.1 


( — ) 
CM 


Kir / fi 


34.35±0.78 


1.242± 0.070 


92% 


4.74± 0.09 


25.6 




Kiririr/fi 


34.83±1.16 


1.407± 0.092 


93% 


4.45± 0.13 


41.8 




Kmr°/fi 


33.69±0.96 


1.191± 0.088 


93% 


4.69± 0.10 


72.0 




Kir/e 


33.75±0.89 


1.212± 0.082 


92% 


4.65± 0.10 


68.8 




Kmnr/e 


33.89±1.32 


1.409± 0.105 


93% 


4.21± 0.14 


73.2 


CM 
O 


Kmr°/e 


35.41±1.06 


1.364± 0.085 


92% 


4.71± 0.12 


93.3 


o 

CM 


Kn/ fi 


32.71±1.06 


1.034± 0.107 


93% 


4.81± 0.12 


37.6 




Krnnr/ fi 


33.84±1.52 


1.323± 0.129 


93% 


4.40± 0.17 


45.6 




Kmr°/fi 


35.52±1.24 


1.406± 0.101 


91% 


4.63± 0.14 


42.9 



4.3 Fit Results 

The fit is performed separately for each data set, providing 18 statistically independent determi- 
nations of ^4i(l)|Vrf,| and of . The results are reported in Table 2. In the same Table, the 
branching fraction, as obtained by integrating the differential decay width, is also shown. The 
agreement between the data and the fit function is good. The distribution of the x 2 probabilities 
is uniform. The measured w distributions are shown in Figure 3, where the three distributions for 
the D° decay modes are summed. 

It should be noted that, because the statistical errors are small, the measurements are dominated 
by the systematic errors, which introduce considerable correlations among the samples. For this 
reason, first the systematic error is discussed, and then the average is presented. 

4.4 Systematic Errors and Consistency Checks 

The individual sources of systematic errors are reported in Table 4 and are described in detail 
below. They can be grouped as 

• Factors affecting the overall normalization (see Equation 3) : Each of these factors is varied 
by its uncertainty, and the resulting variation in the parameters is added to the systematic 
error. The uncertainties on Nr, foo, t b o and B{D* + ) are completely correlated for all the 
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Figure 3: The experimental distribution of w compared to fit results separately for electrons (left) 
and muons (right). The plots from all the different D° decay modes and years of data taking are 
summed together. The points are the data and the histograms are the results of the fit to the sum 
of the signal and the different sources of background. The fit residuals are also shown below. 



samples. The uncertainties on the D branching fractions depend on the decay mode and are 
only partially correlated [18]. 

• Experimental uncertainties in the efficiency: These include track reconstruction, lepton and 
kaon identification, 7r° selection, vertex reconstruction and the cut on the Dalitz plot for the 
K~it + it mode. 

For high-momentum tracks, the tracking efficiency is determined comparing the independent 
information from the SVT and DCH. This approach results in ±0.8% systematic uncertainty 
in the reconstruction efficiency for each track; the errors from each track are added linearly 
(i.e., the error on the branching fraction is N times the single track error, where N is the total 
number of tracks employed, including the lepton). Most of the vr s i ow do not reach the DCH, 
and are reconstructed in the SVT alone. The efficiency for these low-momentum tracks is 
computed from the angular distribution of the vr s i ow in the D* + rest frame. To evaluate this, a 
large set of D* + — ► , D° — > K~ir + decays is selected from generic hadronic events. For 

fixed values of the D* + momentum, the observed angular distribution is compared to the one 
predicted for the decay of a vector meson to two pseudoscalar mesons, and the deviation from 
the expected shape is attributed to the inefficiency of track reconstruction. This inefficiency 
is parametrized as a function of the particle momentum in the laboratory frame. The values 
of the parameters describing the efficiency are computed in different bins of the polar angle 
of the track, both in the data and in the simulation. The systematic error is computed by 
varying each parameter of the efficiency function by its uncertainty, including correlations. 

The efficiencies and systematic errors for lepton and kaon identification and 7r° reconstruction 
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are measured from control samples. The error on lepton identification is common to all 
samples (separately for each lepton type). Looser kaon identification criteria are imposed on 
K~ir + decays; the corresponding systematic error is therefore smaller than for K~ir + ir + ir~ 
and K~-k + tt decays. 

The systematic errors due to the cuts on the vertex probability and on the Dalitz plot are 
evaluated by varying these cuts. 

The efficiency corrections in data and in Monte Carlo simulation are largely independent of 
w, and therefore the value of p^ is scarcely affected, while a substantial error is induced 
on „4,i(l)|Vd,| and branching fraction. Not surprisingly the most noticeable exception is the 
7r s i ow tracking efficiency correction, which affects both *4i(l)|Vd>| and p\ . 



Table 3: Decay modes used in the Monte Carlo simulation of the B — > D* + £~ v^X . These branching 
fraction values are based mainly on theoretical estimates. The theoretical models adopted to 
generate these events are also listed. For the four-body B — > D* + £~i7£ir decays, the Goity and 
Roberts model (GR) [19] was used, while the Isgur and Scora model (ISGW2) [20] was adopted for 
the resonant B — > D**£ + U£ decays. 



B decay mode 



B{B) 
xlO" 2 



Model 



D** decay 



B(D) 



Overall BR 
xl(T 4 



B° 

B°- 

B°- 

B° 

B°- 



D*-Tr u £ + u e 
-> D^£+V(, 
Dl-£ + is e 
+ D*~£ + u e 
-> D'-£ + u e 
* D*'-£+v e 



0.10 
0.56 
0.37 
0.37 
0.02 
0.22 



GR 
ISGW2 
ISGW2 
ISGW2 
ISGW2 
ISGW2 





-> D*~ 


7T° 


0.33 


D*~ 


-► D*- 


-7T° 


0.33 


D* 2 - 


-> D* 


-7T° 


0.103 


D'- 


-> D*- 


-7T 


0.33 


D*'- 


-> D* 


-7T 


0.17 



10. 
18.5 
12.2 
3.81 
0.67 
3.74 



B+ — ! 
B+ 
B+ 
B+ 
B+ 
B+ - 



D*-TT + £+u e 

-»• DJ ^^ 
-»• £>|°^ + ^ 

D*'°£ + v e 



0.20 
0.56 
0.37 
0.37 
0.02 
0.22 



GR 
ISGW2 
ISGW2 
ISGW2 
ISGW2 
ISGW2 







^+ 


0.67 


Df 


-» D*~ 


-vr+ 


0.67 


Df 


-> D*" 


-vr+ 


0.21 


jyo 


-> Z?*" 


-7T+ 


0.67 




-> D* 


-vr+ 


0.33 



20. 
37.5 
24.8 
7.78 
1.32 
7.26 



• Background subtraction: The fraction of background events is determined in data for each 
w bin as described above. This procedure considerably reduces the systematic error. The 
statistical error on the background is automatically accounted for by the fit (see Equation 7). 
A residual model dependence is considered for the D** background. Including narrow, wide 
and non-resonant D* + ir states, a total of twelve modes (for B° and B + semileptonic decays) 
is considered. The D** background is computed assuming the branching fractions reported in 
Table 3. In order to estimate the systematic error due to the uncertainty in the shape of the 
D** background, the cos0Bo t D*e fits m the different w bins are performed considering only 
one decay mode a time, and the corresponding values of \V c b\ and of p^are then determined 
again. The systematic error is then computed as half the difference between the maximum 
and minimum values of the parameters obtained from this study. 

To test whether the background from uncorrelated events is properly handled, the cut on 
cos 6d*£ < is removed and the study is repeated. This results in an increase of the uncorre- 
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lated background by a factor of about 2 in the full cosO B o D *£ range, and by 30% in the signal 
region. No appreciable change in the fit result is observed. As an alternative, the fraction of 
uncorrelated events is determined by counting the events in the rejected region cosOdh > 
and propagating this number to the complementary region using the Monte Carlo. The frac- 
tion obtained is then fixed in the cos6 B o D *£ fit, which is used to determine the fraction of 
D** and signal events only. Once again, no appreciable change in the result is observed. 
Finally, the hypothesis that the ratio of the amount of D** and uncorrelated background over 
the signal be the same for all D* + decay modes is checked with the simulation. While the 
test is satisfactory for the uncorrelated background, a slight discrepancy is observed for D** 
events. For K'ir + 7r + 7r~ decay the fitted fraction is higher by 30% than for the other modes. 
The effect is the same for electron and muon events and does not depend on w. Therefore 
the w fit is repeated by increasing the D** amount in the K~tt +/ k +/ k~ sample and decreasing 
it correspondingly in the K~ir + and K~tt + tt . The difference in the result is propagated as 
systematic error. 

• Theoretical uncertainties: A considerable fraction of the allowed phase-space is removed by 
the lower limit on the lepton momentum. This induces a model-dependence in the computa- 
tion of the branching fraction. The lepton spectrum is determined partly by the fitted shape 
of the Isgur-Wise function, and partly by the angular profile of the decay, which is governed 
by the parameters R± and R2 defined above. Even with perfect acceptance, the uncertainty 
on Ri and R2 directly affects p\ and «4i(l)|Kb|) because these parameters enter in the fit 
function. The values of these parameters determined by the CLEO collaboration [21] are 
used in this analysis. The systematic error is taken as the observed variation of „4i(l)|V^,|, 
p\ and of B(B° — > D* + £~i> i ), when R± and R2 are floated by their uncertainties, accounting 
for their correlation. 

• Fit method: The fit procedure is validated using an independent sample of Monte Carlo 
simulated events, of the same size as the data, and following the procedure described above, 
using however a linear extrapolation for the Isgur-Wise function. The fitted values of p\ x and 
«4i(l)|Vcb| are consistent within their statistical errors with the input values. The statistical 
error on this test is added to the systematic error. In addition, to test that the fit method does 
not bias the result, a set of 500 toy experiments are performed, in which random events are 
generated with realistic efficiency and resolution. Each toy experiment has the same sample 
size for data and Monte Carlo as in the actual measurement. Within errors, no difference is 
observed between the average of the fitted parameters and their generated values. Also, it is 
observed that the pull width is consistent with one. 

4.5 Results 

The results from all subsamples are combined by using the COMBOS package [22], taking into 
account the correlations between the samples. COMBOS was originally developed by the LEP- 
B-Oscillation working group to combine the results of AM^ from different experiments, and was 
then adapted by the LEP-V c b working group to compute the world average for |V^,|, pj^ and 
B(B° — > D* + t~vi). Besides computing averages, errors, and confidence levels, COMBOS also 
provides a breakdown of the single error sources, which is reported in Table 4. 

All the parameters affecting the normalization are treated as fully correlated. The errors on 
the branching fractions of the D° to each final state are fully correlated for each state, but, at 
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Table 4: Summary of all statistical and systematic errors. 



error contribution 


Mi(i)|y c6 |Mi(i)|y c6 | (%) 




5B/B (%) 


statistics (data and MC) 


0.7 


0.02 


0.8 


particle identification 


0.5 


_ 


0.9 


7r slow efficiency 


1.3 


0.02 


1.9 


tracking & 7r° efficiency 


1.3 


- 


2.7 


D* + i~ vertexing efficiency 


0.5 


- 


1.0 


w fit method 


0.6 


0.02 


1.2 


cos9 B o D *£ x 2 fit binning 


0.5 


- 


1.0 


D** background composition 


1.8 


0.06 


2.0 


total number of B produced 


0.6 


- 


1.1 


I (46 j rest frame B momentum 


0.3 




0.7 


7?i (I) and Bo(l) 


1.8 


0.27 


1.8 


total systematic error 


3.4 


0.28 


5.0 


T B 


0.6 






B(D°) 


1.1 




2.0 


B(D*+ D°tt + ) 


0.4 




0.7 


B(T(4S) 


1.3 




2.7 


total systematic error 


1.8 




3.5 



present, the correlations among different decay modes are neglected. The systematic errors on 
particle identification are considered as fully correlated among data taking periods. The error 
on tracking is determined by systematic effects common to all years and therefore is considered 
as fully correlated among all the samples. In contrast the error on 7r s i ow tracking efficiency is 
still dominated by the statistical uncertainty from the control samples; therefore it is completely 
correlated for samples collected in the same year of data acquisition, but uncorrelated for different 
years. Model errors are common to all data sets. 

The mean values are obtained by first averaging results for different running periods, and then 
computing the average result for the three years. The confidence level of the «4i(l)|Vd,| and p\ x 
average values is 30%. The results are reported in Table 5. 



Table 5: Fit results for the three running periods. 



year 


A^Wcl x 10 3 p% 


2000 
2001 
2002 


33.60±0.46±1.33 1.23±0.04±0.28 
34.63±0.36±1.44 1.25±0.03±0.28 
34.73±0.45±1.39 1.31±0.04±0.28 


Total Average 


34.03±0.24±1.31 1.23±0.02±0.28 



For the sum of all 18 data subsamples the following result is obtained: 

^i(l)|F c6 | = (34.03 ± 0.24 ± 1.31) x 10" 3 , (9) 
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p\ = 1.23 ±0.02 ±0.28, 
B(B° D* + £-i7 e ) = (4.68 ± 0.03 ± 0.29) %, 



(10) 
(11) 



where the first error is statistical and the second is systematic. The statistical correlation between 
^li(l)|y c6 |and p^is 92%. 

5 Conclusions 

A sample of 55,700 signal events from the decay process B° — > D* + l~vi is selected from a set of 
86 x 10 6 T(AS) -> BB decays collected by the BABAR detector. The product of \V cb \ and the 
form factor at zero recoil, = A\{\), is measured, as well as the derivative of the form factor, 
p\ v again at zero recoil. The integrated branching fraction is also computed. These results are in 
agreement with those obtained by Belle [8] and with the LEP average [18], but differ by more than 
two standard deviations from the CLEO result [9]. Using the value -4(1) = 0.913lo!o35 reported 
in [4], based on a lattice QCD calculation, the value: 

\V ch \ = (37.27 ± 0.26(stat.) ± 1.43(syst.)+};23(th.)) x 10" 3 (12) 

is then computed. 

From the measurements of the hadronic mass moments and decay rate for the inclusive semilep- 
tonic B decays BAEL4R has obtained another determination of \ V c b\ [23]. The preliminary result is: 

\V cb \ = (42.1 ± 1.04(exp.) ±0.72(th.)) x 10~ 3 (13) 

Both the experimental and theoretical errors are independent. By interpreting the theoretical errors 
as ±68% confidence interval, the two values differ by 2.0 standard deviations. 
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